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Abstract 

We consider non-relativistic electrons in one dimension with in- 
finitely strong repulsive delta function interaction. We calculate 
the long-time, large-distance asymptotics of field-field correlators 
in the gas phase. The gas phase at low temperatures is character- 
ized by the ideal gas law. We calculate the exponential decay, the 
power law corrections and the constant factor of the asymptotics. 
Our results are valid at any temperature. They simplify at low 
temperatures, where they are easily recognized as products of free 
fermionic correlation functions with corrections arising due to the 
interaction. 
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Introduction 

Non-relativistic electrons in one dimension, interacting strictly locally through 
a delta potential constitute one of the important exactly solvable many-body 
systems with applications in solid state physics. Historically this was the 
first system to be solved exactly by nested Bethe ansatz |], |J. Soon after an 
exact formulation of its thermodynamics was obtained B. In this article we 
address the problem of the exact calculation of its correlation functions. Our 
work is based on a recently derived determinant representation for the field- 
field correlators [f|, || valid in the case of infinite repulsion. The determinant 
representation gives the field-field correlators as functions of space and time 
variables x and t in a grand canonical description, i.e. the correlators depend 
on temperature T, chemical potential /i and on the external magnetic field 
B. We study the the long-time, large-distance asymptotics of the field-field 
correlators. We shall consider the limit t — > oo for fixed ratio k = x/2t. 

Let us write the Hamiltonian in terms of canonical anticommuting quan- 
tum fields ip a (x, t) (a =|, |), 

H = |" dx {a^^Va + c : (ipl^af : - /x ViV'a + B(VVV)} • (1) 

Here a z is a Pauli matrix, and c is the coupling constant. In the following we 
shall consider infinite repulsion c = +00. Then there is an explicit expression 
for the pressure as a function of temperature, chemical potential and magnetic 
field §, 

P = — dk\n(l + e ^ +B - k ^ T + e ^ B ~ k ^ T ) , (2) 

which may serve as thermodynamic potential. Note that this expression is 
formally the same as for a free Fermi gas with effective chemical potential 
f3T = fi + Tln(2ch(5/T)). Hence there are two different zero temperature 
phases depending on whether lim T ^ 0+ (3T = p,+ \B\ is positive or negative. If 
[A + \B\ > 0, the density D = dP/dfi has a positive limit as the temperature 
goes to zero. If /i + \B\ < 0, the density at zero temperature vanishes. This 
is the case we are interested in. We shall assume throughout the remainder 
of this article that (3 < 0. We call the phase determined by this condition the 
gas phase, since in this phase pressure and density at low temperatures are 
related by the ideal gas law, 



P = DT 



(3) 
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As the temperature increases the density becomes large, and the ideal gas 
law does not hold anymore. There is a macroscopic number of electrons in 
the system, which may carry a current. At strictly zero temperature there 
are no particles in the system, and the system is an insulator. 

We shall consider the following two-point correlation functions of local 
fields, 

_ tr(e-^(x,^t(o,o)) 
^ a p\ x ^) ~ ti{ e - H l T ) ' 

tr (e- g / r V>t(s,t)^(0,0)) 
(*ap{ x ,t) - tr(e- H / T ) ' ^' 

where a, (3 =|, |. Due to conservation of the total spin G+Ax, t) = GfAx, t) = 
0. Furthermore, due to the invariance of the Hamiltonian (p]) under the trans- 
formation B# —B, we have the identity G^(x,t) = G^(x,t)\B^-B- 
Hence, we can restrict our attention to the correlation functions G*^(x,t). 



Determinant representation 

Let us recall the determinant representation for the correlation functions 
G++(x, t), which was derived in [||, |5|. We shall basically follow the account of 
|J. Yet, it turns out to be very useful for further calculations to rescale the 
variables and the correlation functions. The rescaling 

x r = -VTx/2 , t r = Tt/2 , (6) 
g ± = G±/Vf , (7) 
h = B/T (8) 

removes the explicit temperature dependence from all expressions. Further- 
more, it will allow us to make close contact with results which were ob- 
tained for the impenetrable Bose gas [7L [5], [|. The index V in (|B|) stands for 
'rescaled'. For the sake of simplicity we shall suppress this index in the follow- 
ing three sections. We shall come back to physical space and time variables 
only in the last section, where we consider the low temperature limit. 



4 



The rescaled correlation functions g + and g in the rescaled variables can 
be expressed as 0, 



_ 2ii(/?-/i-ln(2ch(/i))) fn „ i ^ ^ 

S + (*,t) = / rfn l7 '7 & ++ det(/ + y) , (9) 

Z7T J-tt 1 — cos(r]) v ' 

-2it( i 9-/l-ln(2ch(/i))) /*7r A 

ff "(a;,t) = — J_J V F(y,r))B^det(I + V) . (10) 

Here 7 and -F(7, 77) are elementary functions, 

7 = l + e 2h , (11) 
e ir) e~ iv 

F(j, V ) = 1 + - + — . (12) 

7 — e iri 7 — e 177 

det(J + V") is the Fredholm determinant of the integral operator / + V , where 
I is the identity operator, and V is defined by its kernel V(\, //). A and /i are 
complex variables, and the path of integration is the real axis. In order to 
define V(\,n) we have to introduce certain auxiliary functions. Let us define 

r(A) = i(A 2 t + Ax) , (13) 

= ' (14) 

poo £-2t(m) 

E{\) = p.v. / dfi- , (15) 

J-00 7l{fl — A) 



e.(A) = , (16) 



7T 



e+(A) = ^^e-^){(l-cos(r ? ))e-W 



1 - cos(r]))e 2rW E(X) + sin(^)} . (17) 

Note that $(A) is the Fermi weight. V(A,/x) can be expressed in terms of e + 
and e_, 

nA^)^ e+(A)e - ( ^~ e+Me - (A) • (18) 

A — yU 

Denote the resolvent of V by -R, 

{i + v)[i-R) = (l-R)(l + V)=i . (19) 
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Then R is an integral operator with symmetric kernel, 

R(x ^ } = UWf-M-fMf-W , (20) 

A — /i 

which is of the same form as V(X,fi). The functions f± are obtained as the 
solutions of the integral equations 

/oo 
d/iV(A,/i)/ ± Oi) = e±(A) . (21) 
-oo 

We may now define the potentials 

/■OO /"OO 

B ab = dXe a (X)f b (X) , C ab = dXXe a (X)f b (X) (22) 



for a,b = ±. B enters the definition of g (x,t), equation (PD|). b ++ in (Q) 

is defined as 

b++ = B ++ -G(x,t) , (23) 

where 

G(x, t) = 1 ~ ^ f dA e~-W = (1 - C ° S ^ e " W4 e- 2 /^ . (24) 

27T J-cxD 2V27Tt 

The remaining potentials i? a b and C a b will be needed later. 

It is instructive to compare the determinant representation (|9|) for the cor- 
relation function g + (x,t) with the corresponding expression for impenetrable 
Bosons (cf e.g. page 345 of ||). The main formal differences are the occurrence 
of the rj- integral in (||) and the occurrence of r] in the definition of e + . As can 
be seen from the derivation of (pj) in ||, the i]- integration is related to the spin 
degrees of freedom. For r] = ±7r the expression — |e 2l/3t 6 ++ det(7 + V) agrees 
with the field-field correlator for impenetrable Bosons (remember, however, 
the different physical meaning of f3). 

Differential equations 

As in case of impenetrable Bosons [0, |J it is possible to derive a set of 
integrable nonlinear partial differential equations for the potentials b ++ and 

B and to express the logarithmic derivatives of the Fredholm determinant 

det(J + V) in terms of the potentials B ab and C ab . 
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The functions f± satisfy linear differential equations with respect to the 
variables x, t, and (3, 

i C:)=*®=*®= ■ (25) 

The Lax operators L, M and N are given as 

L = d x + iXa z -2iQ , (26) 

M = d t + i\ 2 a z - 2iXQ + d x U , (27) 

N = 2X8,3 + d x + 2\tXo z + ixa z - AitQ - 2d p U , (28) 

where the matrices Q and U are defined according to 



I o b ++ | 














\b 






B+- 



Q=\ , U=\ T ^ . (29) 

Mutual compatibility of the linear differential equations fl25|) leads to a 

set of nonlinear partial differential equations for the potentials b ++ and B 

Instead of using b ++ and B , it is convenient to introduce the following 

notations (cf. p. 352 of the book ||), 

g+ = e^b ++ , g_ = <r*# B _ , (30) 

n = 9+9- = K+B— , p = g-d x g + - g+d x g- . (31) 

The space and time evolution is driven by the separated nonlinear Schrodinger 
equation, 

-id t g + = 2Pg + + ld 2 x g + + 4g 2 + g„ , (32) 
idtg- = 2Pg_ + ±d 2 x g„ + 4glg + . (33) 
p and n satisfy the continuity equation 

-2\d t n = d x p . (34) 
The equations containing derivatives with respect to j3 are 

-idt<p + 4df,p = , (35) 
d x (p + 8dpn + 2 = . (36) 
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Here the function ip is defined as 

/ . „s dpd x g+ dpd x g- 

<P{x,t,P) = = . 37 

9+ 9- 

To describe the correlation functions (§) and fllPf ) one has to relate the 
Fredholm determinant det(/ + V) and the potentials B a b and C a b. Let us use 
the abbreviation a(x,t,/3) = In det (/ + V") . The logarithmic derivatives of the 
Fredholm determinant with respect to x, t and (3 are 

d x a = -2LB+_ , (38) 
d t a = -2i(C+_ + C_ + + G(x,t)B__) , (39) 
dpa = -2i%(C+_ + C_+ + G(x, t)B__ ) - 2ix9 /JJ B + _ - 2(<9 /3 B+_) 2 

-2it(B__^6 ++ -6 ++ a^__) + 2(^6 ++ )(^B__) . (40) 

For the calculation of the asymptotics of the Fredholm determinant we further 
need the second derivatives of a with respect to space and time, 

d 2 x a = AB__b ++ , (41) 
d x d t a = 2i( J B__Q B & ++ - b++d x B—) , (42) 
d\a = 2i(B^d t b ++ -b ++ d t B__) + 8B 2 _b 2 ++ + 2(d x B^)(d x b ++ ).(43) 

Note that 

lim a = . (44) 

This follows from linig^-oo i?(A) = and is important for fixing the integration 
constant in the calculation of the asymptotics of the determinant. 

The remarkable result of this section is that the differential equations for 

b ++ and B as well as the expressions (|38|)-(l40|) for the logarithmic derivative 

of the Fredholm determinant are of the same form as in case of impenetrable 
Bosons 0,0. 

Asymptotics of the correlation functions 



The differential equations (p2|), fl33|) have many solutions. b ++ and B are 

the particular solution of (|32]) and (j33|) which is fixed by the integral equations 
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(|21"p. Alternatively b ++ and B can be obtained by solving a Riemann- 

Hilbert problem, which was formulated in 0. The Riemann-Hilbert problem 

is the appropriate tool to obtain the leading asymptotics of b ++ , B , d x a and 

d t <y. In fact, using the method first suggested by Manakov [|K| and further 
developed in and in ||, it was shown in that the following estimates 
hold, 

6 ++ = 0(H) , £L_ = fl(H) , (45) 

and 

d x a = - r d\sign(\-\ )Htp(\,fl)) + O(r 1 *) , (46) 
d t a = - r dX Asign(A-A )ln(^(A,/?)) + e>(H) . (47) 

The function ip(\, (3) in the above expressions is defined as 

(p(\, (3) = 1 + 0(A) ^-^ignix-xo) _ ^ _ ( 4 g) 

A is the stationary point of the phase r(A), t'(A ) = 0, i.e. 

Ao = -| • (49) 

Equations (|45| ) - (|49|) came out of the direct asymptotic analysis of the 
Riemann-Hilbert problem. As we shall explain in the following these equa- 
tions, together with fl38l)-(fi~3|), fix the complete asymptotic expansion of the 
correlation functions g + (x,t) and g~(x,t) up to a single unknown function. 

It follows from (|45[) that b ++ and B are a decaying solution of the sep- 
arated nonlinear Schrodinger equation fl32f), fl33[ ). The form of the complete 
asymptotic decomposition of the decaying solutions of the separated nonlinear 
Schrodinger equation is known (12], |13|| , 

j / - 2n ln fc 4t \ fix 2 . , 1 

K+ = t 2 ^u + ^ —^r u nkj exp j— - izy In 4t| , (50) 

5— = ^ + ^^__ Wnjfc Jexp|-— +i^ln4t| . (51) 

Here the quantities «o, t>o, «nit, t> n fc and z/ are functions of Ao = —x/2t and 
of /3 and rj. Inserting the above asymptotic expansions for B and b ++ into 
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the differential equations fl32|), (|33|) we obtain expressions for u n k, v n k and v 
in terms of Mo and Vq. In other words, the two unknown functions Mo and t>o 
determine the whole asymptotic expansion (]5Ti|), (|5"T|). In particular, 

v 

UoV = -- ■ (52) 

This relation can be used to calculate v. fl4"T|), (|50D, (|51|) and ( |52"D imply that 

n , z/ ,/ln 2 4t\ . , 

d 2 x a = — t +0 l^—j . (53) 

Differentiating (|4*6| ) with respect to x and comparing the result with the right 
hand side of (|53f) we find 



v = ~ In (1 - 2(1 - cos(t 7 ))i?(Ao)(1 - 0(A„))) , (54) 

such that the product u v is fixed, and we are indeed left with only one 
unknown function, say, u . 

The second logarithmic derivatives of the Fredholm determinant with re- 
spect to space and time are given as functions of b ++ and B by the equa- 
tions d4lD-(^3|). Hence, their complete asymptotics, too, is determined by 
the expansions ( |50D and (|5T|). We may integrate (|4"I|)-(fl3"|) to obtain d x a and 



dto. The integration constant, which is a function of /3, is fixed by the leading 
asymptotics (pjj), (f|7|). Then, integrating fl3"5|)-(|4"U[) yields a up to a numerical 



constant, which follows from the asymptotic condition (f4|). The calculation 
is the same as for the impenetrable Bose gas and can be found on pages 455 
- 457 of the book ||. It results in the asymptotic formula 

1 r°° v 2 v 2 

a = - dX\x+2Xt\\n((p(X,(3)) + — hi4i+— - 

7T J-co 2 2 

+ 2i / d(3 (u dpv - v dpu Q ) 

J —CO 

+ i lljP { dp /" dAsi S n ( A - A o)ln(^(A,/3)) 
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So far the equations (pD|), (p|), (|52"D, ( |53"D and ( |55D determine the asymp- 
totics of the correlation functions g + (x,t) and g~(x,t) up to an unknown 
function Uq. In order to calculate u$ one has to reexamine the Riemann- 
Hilbert problem. The asymptotic calculation of || has to be refined. For 
impenetrable Bosons a similar refinement was achieved in ||. Fortunately, 
the result of || depends only on some general features of the functions enter- 
ing the Riemann-Hilbert problem in its canonical form. It turns out that it 
also applies in our case. We can use equation (6.21) of M to obtain 



u ° = V 2 2^(A )r(-i,) (56) 

= i|sin(r ? /2)| v ^exp((A^-/3)/2 + i^ ) , (57) 



where 



3tt 

= -— + ar g r(iz/) + , (58) 
= -- r c /Asign(A-A )ln|A-A |9 A ln(^(A,/?)) . (59) 

7T J-oo 

An alternative rigorous derivation of these formulae can be performed via the 



Deift-Zhou nonlinear steepest descent method |L4j . It is also worth comparing 
equations fl57D-(|59D with the asymptotic formula for the solution of the Cauchy 
problem of the NLS equation obtained by Zakharov and Manakov (see 



also 021, 131 and [16 



Let us summarize our results. Inserting fl5"0|), ([ST]) , (|55| ) and (^7[) into (|9|) 
and (pTDj) we obtain the following expressions for the leading asymptotics of 
the correlation function, 

g+( X ,t) = e il ' 2 /2i+2it/3 e -2it(ft+ln(2ch(h))) 

r ^1L C + (a,, A „>(«)*<->' 

1 — cos(r/) 

exp(- /°°c/A|x + 2At|ln(^(A,/?))l , (60) 

I 7T J-oo J 
g-( X ,t) = e -^ 2 /2<-2it/3 e 2ii(/ i +ln(2ch(/ l ))) 
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d V ^^-C-(X ,M (4t)^ +i > 2 

27 



exp(- H d\\x + 2\t\ln{ip{\,(3))\ , (61) 

I 71 J -oo J 



where 



v 2 



C + (X ,(3 }V ) = -|sin(r 7 /2)|^-exp|-(A^-/3)+i* + y 



f d(3 (xp/2 + vdpVo) (62) 
^/ (fy y_ oo rfAsign(A-A )ln(^(A,/3)) 



2tt 2 

C-(A 0) /3,r/) = C + (Ao,/3,r ? )exp(-(A2-/3)-2iM/ )/sin 2 (r ? /2) . (63) 

These equations are valid for large t and fixed finite ratio Ao = —x/2t. Cor- 
relations in the pure space direction t = will be discussed in a separate 
publication. We would like to emphasize that (|60D and ([H]) hold for arbi- 
trary temperatures. Note that there is no pole of the integrand at rj = 0, 
since ^fv ~ \r]\ for small r\ and thus C + {\q,(3, rj) ~ rf. 



Asymptotics in the low temperature limit 

For the following steepest descent calculation we transform the //-integrals 
in (|60D , (|BTD into complex contour integrals over the the unit circle, setting 
z = e ir? . Since we would like to consider low temperatures, we have to restore 
the explicit temperature dependence by scaling back to the physical space 
and time variables x and t and to the physical correlation functions Gp*. 
Recall that in the previous sections we have suppressed an index V referring 
to 'rescaled'. Let us restore this index in order to define k = XqVT = 
x/2t, &{k) = Mk/VT), <p(k,P) = cp r {k/VT,[3), C±{k ,/3,z) = C±{\ ,P,r)), 
F(j,z) = F r (~t,ri). Then 

G+{x,t) = 2iVTe ix2/u+it ^- B) 

idz F( - J '*l c + (k , P, z)(2?Y)£W»)-i) 3 e *gW , (64) 
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G n (x,t) = -iVT e -^-^-m 



dz^^C-(k ,(3,z)(2Tt)^ z ^ 2 e ts ^ , (65) 
272; 

1 r°° 

S( z ) = - dk\k-k \ln(cp(k,P)) . (66) 

71 J-00 



where 



We would like to calculate the contour integrals (|64|) , fl65|) by the method of 
steepest descent. For this purpose we have to consider the analytic properties 
of the integrands. Let us assume that k Q > 0, and let us cut the complex 
plane along the real axis from —00 to — e _/3 and from —e@~~ k o/ T to 0. The 
integrands in (|64]) and ( |65|) can be analytically continued as functions of z 
into the cut plane with the only exception of the two simple poles of -F(7, z) 
at z — 7 ±1 . We may therefore deform the contour of integration as long as 
we never cross the cuts and take into account the pole contributions, if we 
cross z = 7 or z = 7 _1 . 

The saddle point equation dS/ dz = can be represented in the form 

dk k poo d/lx, \z 



1 _|_ ze -P e {k+k y/T 



(67) 



This equation was discussed in the appendix of 0. In it was shown that 
( |B"TD has exactly one real positive solution which is located in the interval [0, 1]. 
It was argued that this solution gives the leading saddle point contribution 
to (l&l) and fl65|). At small temperatures (|67D can be solved explicitly. There 
are two solutions z± = ±z c , where 

e- fc o/2T (68) 



C 2n^k 3 /2 

In the derivation of fl68|) we assumed that ko 7^ 0. The case ko = has to be 
treated separately (see below). 

The phase tS(z) has the low temperature approximation 

tS(z) = -2k Dt I (l - zl + 1 - z\ . (69) 

Here D = dP/d/i is the density of the electron gas. The low temperature 
expansion (E3) is valid in an annulus e^ _fc o/ T <^ \z\ <^ e _/3 , which lies in our 
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cut plane. The unit circle and the circle \z\ = z c are inside this annulus. We 
may thus first apply (|69D and then deform the contour of integration from 
the unit circle to the small circle \z\ = z c . Let us parameterize the small 
circle as z = z c e ia , a G [— ir, 7r]. Then S(z(a)) = —2k D((z c — l) 2 + 2z c (l — 
cos(a))), which implies that the small circle is a steepest descent contour 
and that on this contour S(z-) < S(z) < S(z+). The maximum of S(z) 
on the steepest descent contour at z = z + is unique and therefore provides 
the leading saddle point contribution to (|64|), (|65| ) as t — > oo. The saddle 
point approximation becomes good when tS(z(a)) = —2k Dt((z c — l) 2 + 
z c a 2 + 0(a 4 )) becomes sharply peaked around a = 0. Hence, the relevant 
parameter for the calculation of the asymptotics of G++ is 2k$Dt = xD rather 
than t. xD has to be large compared to z~ l . The parameter xD has a simple 
interpretation. It is the average number of particles in the interval [0, x] . Let 
us consider two different limiting cases. 

(i) xD — > 0, the number of electrons in the interval [0, x] vanishes. In 
this regime the interaction of the electrons is negligible. An electron 
propagates freely from to x. G++ cannot be calculated by the method 
of steepest descent. We have to use the integral representation (|60D , (|61| ) 
instead. Since tS(z) and u(z) tend to zero on the contour of integration, 
the integrals in (|60D and ( |6l]) are easily calculated. We find 

GftCM) = 0r^^ B »e^/ 4! , (70) 

As expected, equations (ffO| ) and (ffll) are the well known result for free 
fermions. 

(ii) xD ^> z' 1 , the average number of electrons in the interval [0, x\ is large. 
This is the true asymptotic region, x — > oo. In this region the interaction 
becomes important. At the same time the method of steepest descent 
can be used to calculate This case will be studied below. 

In the process of deformation of the contour from the unit circle to the 
small circle of radius z c we may cross the pole of the function Ffa, z) at 
z = 7 _1 . Then we obtain a contribution of the pole to the asymptotics 
of Gf+. It turns out that the pole contributes to G^, when the magnetic 
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field is below a critical positive value, B c = ky/4. Below this value the pole 
contribution always dominates the contribution of the saddle point. Hence, 
we have to distinguish two different asymptotic regions, B > B c and B < B c . 
On the other hand, if we consider the asymptotics for fixed magnetic field, 
we have to treat the cases B > and B < separately. For B > we have 
to distinguish between a time like regime (k^ < AB) and a space like regime 
(/cq > 4B). In these respective regimes we obtain the asymptotics: 

Time like regime (x < 4t\Z~B): 

Gt(x,t) = A . 1 r i-M*c) e i*fr-B) isV4t -sg, , 72) 

G7Jx,t) = A ^; B - k ° )/T r ^) e -«M) e -^/^ e -*Di (73) 

where 

= - ^574 ' ^ 

Di = iffiW . (75) 

Di = dPj d(fi + B) is the low temperature expression for the density of down- 
spin electrons. 



Space like regime (x > AtyB): 

G+(x,t) = r \^ {l ^)^-B) e ^/U e - X D^ (?6) 
G^(x,t) = ^ e (^-feo 2 )/^ r l+M7- 1 ) e -i^- J B) e -- 2 /4t e -^ i)(77) 



where 



( 3B+M _fc2 )/T 

K7- 1 ) = "— -= ■ (75 

Z7T 



For B < there is no distinction between time and space like regimes. The 
asymptotics is always given by the equations (|76| ) and ([77]). 

It is instructive to compare the asymptotic expressions (|7^)-fl7Tj|) for the 
correlation functions with the corresponding asymptotics for free fermions, 
fl7g), flnp. We can interpret the factors r^^e" 2 ^ /JAttz c xD 1 in (J72|) and 
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( |73|) and t ±w(cl 1 ^e~ xD - 1 in (|76| ) and ([TT]) as low temperature corrections to the 
free fermionic correlation functions. The occurrence of the density of down- 
spin electrons in the exponential factors in equations (]72"|)-(|Tg|) has a natural 
interpretation. Correlations decay due to interaction. Because of the Pauli 
principle and the locality of the interaction in our specific model, equation 
([!]), up-spin electrons only interact with down-spin electrons. Therefore the 
correlation length is expected to be a decreasing function of the density of 
down-spin electrons, which diverges as the density of down spin electrons goes 
to zero. We see from (|72|)-(|78D that the low temperature expression for the 
correlation length is just and thus meets these expectations. 

In the limit B — > — oo, fi — > — oo, fi — B fixed there are no J,-spin electrons 
left in the system, — > 0. This is the free fermion limit. In the free 
fermion limit B < B c , and the asymptotics of G^(x,t) and Gj^(x,t) are 
given by equations (|7E|), ([77D, which turn into the expressions ([7[]), Q7T| ) for 
free fermions. 

The pure time direction ko = requires a separate calculation. For k$ = 
the saddle point equation (p7|) has the solutions z = ±1 for all temperatures. 
The unit circle is a steepest descent contour with unique maximum of S(z) 
at z — 1, which gives the leading asymptotic contribution to the integrals in 
(|64D and (|65|). We find algebraically decaying correlations, 



G+(0,t) 



c ,+ r i e tt0x-B) 



G n (0,t) 



(79) 



where 



'l + 2e- 2B ' T ) 



e Qi+B)/T + g(M~ B)/T\ + e Qt+B)/T + e (M-B)/T^ 



(80) 



Cn 



e if i + 2e- 2B ' T 



e (n+B)/T + e <jJ,-B)/T 
1 4. e (v+B)/T + e (^~B)/T 



These formulae are valid at any temperature. 
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